Self-difFusion in binary blends of cyclic and linear polymers 

Sachin Shanbhag^'^'Q 

^School of Computational Science, Florida State University, Tallahassee, FL, 32306 
^Department of Chemical and Biomedical Engineering, 
FAMU-FSU College of Engineering, Tallahassee, FL 32310 
(Dated: May 22, 2008) 

Abstract 

A lattice model is used to estimate the self-diffusivity of entangled cyclic and linear polymers in 
blends of varying compositions. To interpret simulation results, we suggest a minimal constraint 
release model for the motion of a cyclic polymer infiltrated by neighboring linear chains. Both, 
the simulation, and recently reported experimental data on entangled DNA solutions support the 
simple model over a wide range of blend compositions, concentrations, and molecular weights. 

PACS numbers: 02.70.Uu (Applications of Monte Carlo Methods), 83.80.Tc (Polymer blends (rheology)) 
Keywords: rings, cyclic, entanglement, blend, diffusion, polymer melt 



INTRODUCTION 



According to classical theories of polymer physics, flexible chains in solution assume 
coiled conformations. As the polymer concentration, or contour length is increased, these 
coils overlap and produce entanglement effects, which include a pronounced retardation of 
molecular mobility. For linear polymers (LP), this transition from unentangled to entangled 
dynamics is marked by a change in the zero-shear viscosity r/o from rj^ M to tjq ^ M^'^, 
and self-diffusivity from ~ M^^ to ~ M"^'"^, where M is the molecular weight. 
Molecular topology has a strong effect on the dynamics of polymers in the entangled state, 
and constitutes a subject of fundamental and industrial interest in polymer physics and 
rheology. Although important gaps in understanding persist, the behavior of entangled linear 
and branched polymers, such as stars, are relatively well-described at small deformation 
rates, using the tube ansatz. Here, the molecular dynamics of a tracer chain enmeshed in a 
matrix of other molecules, are formulated in terms of a diffusion problem that describes the 
motion of chain ends in a hypothetical tube 

Concentrated solutions of ring or cyclic polymers (CP), which lack chain ends, are scien- 
tifically intriguing, since they defy a simple description in terms of the tube model. While 
interest in CPs has recently been rekindled, the conformational^and dynamic properties 
of CPs in gels and in melts have been studied theoretica 
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ly ^, ^, y, [7 1, computationally 



19[ since the 1980s. Almost all 



the theoretical and computational studies have investigated the characteristics of pure CPs 



and LPs, while blends of CPs and LPs have escaped the same level of scrutiny [111 . I20l . |21 |. 
These blend systems are important for two reasons: (i) most experimental data on pure 
CPs are, in fact, data on cyclic-linear blends (CLE), due to contamination or limitations of 
purification methods, and (ii) the dynamics of such CLBs, are extraordinarily sensitive to 



the concentration of LPs, as demonstrated by the linear viscoe 
CLBs 
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23| . and self-diffusion studies of DNA solutions 
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astic response of polystyrene 



25| . These studies indicate 



a dramatic change in the mobilty of CPs that is, both, unexpected and unexplained, and 
might supply deep insights into entangled polymer dynamics. 

In this letter, a lattice-based dynamic Monte Carlo method called the bond-fluctuation 
model (BFM) is used to monitor the trajectory of LPs and CPs in entangled CLBs. A 
minimal model is constructed, to interpret the diffusivities obtained from the simulation. 



This approximate theory is able to reasonably explain recent experimental data on entangled 
DNA solutions. 



MODEL AND METHODS 

We use Shaffer's version of the BFM 
of CPs and LPs 
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201, 



27l | which has recently been applied to blends 



2l| . Since the model has been described in detail previously, only a 
brief summary presented here. Monomers or beads are placed on a simple cubic lattice. 
To generate an equilibrated CLB, we insert ric non-concatenated CPs and rii LPs, each 
consisting of monomers, on a 3D cubic lattice in a simulation box of size L x L x L. 
To simulate melt-like behavior, the total fractional occupancy of the lattice is maintained 
a.t (j) = (pc + (pL = 0.5, where (pi = riiN/L^ represents the fractional occupancy of CPs 
{i = C) and LPs {i = L). A trial move is attempted by displacing a randomly selected bead, 
belonging to either a CP or LP, by one lattice unit. It is accepted if it does not violate the 
excluded volume, chain connectivity and chain uncrossability constraints. One Monte Carlo 
step (MCS) consists of {til + nc)N trial moves. Throughout this letter, length is expressed 
in units of lattice spacing and time in MCS. 

In this study, we considered two series of CLBs, viz. A^ = 150, and A^ = 300. We varied 
the composition of the blend from cpc = 0.5 to (pc = 0.0, corresponding to the range between 
pure CPs to pure LPs, respectively, as summarized in Table [B The duration of the simu- 
lation, Tsim was picked to ensure that molecules had diffused at least two radii of gyration. 
During this period, we tracked the positions of the molecules by taking snapshots at periodic 
intervals. To determine the self-diffusivity, we calculated the mean-squared displacement of 
the center-of-mass of LPs and CPs separately, via, g^lt) = {{rcm{t + t) — rcmlT))^)- Here, 
the average extends over molecules of a given topology (LP or CP), and r is a dummy 
time variable. The diffusion constant D was calculated from the slope of the mean-squared 
displacement according to the relation, dg^{t)/dt = 6D. 

RESULTS 

To estimate the diffusion constants we performed linear regression analysis on the mean- 
squared displacement gslt) in the interval t = O.lbTsim — O.TTsim, where the function gslt) 



A R t / gsj-Tsim) D / 93( 



N=150 



0.500 


7.95 ± 0.10 












0.479 


8.21 ± 0.07 


3.32 


5.85 


± 


0.16 


4.09 


n /ICO 

0.458 


8.10 ± 0.07 


3.30 


5.82 


1 

± 


0.11 


3.63 


0.438 


8.19 ± 0.08 


3.23 


5.76 


± 


0.10 


3.65 


0.375 


7.99 ± 0.08 


3.15 


5.80 


± 


0.07 


4.37 


0.313 


8.02 ± 0.09 


2.83 


5.50 


± 


0.05 


3.97 


0.250 


8.21 ± 0.10 


2.39 


5.54 


± 


0.05 


3.43 


0.188 


8.16 ± 0.12 


2.32 


5.46 


± 


0.04 


3.90 


0.125 


8.15 ± 0.15 


2.50 


5.30 


± 


0.03 


4.59 


0.063 


8.01 ± 0.18 


2.58 


5.25 


± 


0.03 


5.57 


0.042 


8.04 ± 0.25 


2.68 


5.25 


± 


0.03 


5.81 


0.021 


8.33 ± 0.34 


2.49 


5.16 


± 


0.03 


6.32 


0.000 






5.09 


± 


0.14 




N=300 


0.500 


11.20 ± 0.17 












0.450 


11.32 ± 0.15 


3.00 


8.60 


± 


0.20 


2.64 


0.375 


11.51 ± 0.18 


2.89 


8.71 


± 


0.20 


2.73 


0.250 


12.27 ± 0.25 


2.68 


8.10 


± 


0.11 


3.41 


0.167 


12.02 ± 0.29 


2.50 


7.50 


± 


0.11 


3.77 


0.125 


11.80 ± 0.31 


2.40 


7.32 


± 


0.10 


4.12 


0.100 


12.37 ± 0.33 


2.70 


7.24 


± 


0.08 


4.86 


0.050 


12.35 ± 0.55 


2.75 


7.10 


± 


0.07 


6.27 


0.025 


11.67 ± 0.56 


2.79 


6.96 


± 


0.05 


7.41 


0.000 






7.02 


± 


0.06 





TABLE I: Description of the systems simulated. Simulation box size -Lf,ox = 60, and total density 



0c + </'L = 0.5. The radii of gyration Rl and Rc are reproduced from ref. [20(] for reference. The 
third and fifth columns show that the molecules have diffused at least two times their radius of 
gyration during the simulation. 



was almost perfectly linear. Data on pure CPs and LPs was obtained from literature [id. l26|. 
Figure [1] shows the variation of the diffusivity as a function of the fraction of the LPs for 

= 150 and = 300, respectively. As the linear fraction increases, the diffusivities of 
both, the CPs and LPs decrease, although Dc drops more sharply. This is particularly 
evident for A^ = 300. Further, that decrease is most pronounced at small 0^. 

Experimental data on entangled DNA solutions (fig. 3 from ref. [25]) suggest that LPs 
diffuse more sluggishly in a CP matrix than in a LP matrix, and that the strength of this 
slowdown increases with A^. This is possibly the reason why the decrease in at small 
(pL is more prominent for A^ = 300. Other studies on polystyrene melts (fig. 2 from ref. 
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FIG. 1: The self-diffusivity of the cycUc (triangles) and hnear (circles) molecules for N = 150 (open 
symbols) and = 300 (filled symbols) at different blend compositions. 



17|) suggest that is independent of the composition of the blend, which appears to be 



true over a wide composition range in figure [H The apparent contradiction in these two 
datasets may be partially reconciled through our findings. In all likelihood, a small fraction 
of the supposed polystyrene CPs in the matrix were contaminated with LPs, 171] increasing 
the actual (pi- It is then conceivable that in the window of observation the decrease in D]^ 
at small (j)^ was not captured. 



Minimal Model 



In order to interpret these simulation results, we suggest a minimal model, which ignores 
prefactors and other numerical details. In this scheme, we visualize a CP in a blend, which 
is threaded by Zc surrounding LPs. We recall from previous simulations that the average 
number of entanglements on a CP, Zc, varies according to the linear fraction as Zc{(t)L) / Zl = 
where Zl ^ is the average number of entanglements on a LP and is independent 



of blend composition 



2l|. 



The threaded LPs restrain the mobility of the CP. As some of the LPs venture out, 
others arrive and form entanglements at the same rate, and the equilibrium structure of the 



melt is not disturbed. Consequently, the primitive path of the CP itself undergoes local 
rearrangement as it relaxes by constraint release (CR) Rouse motion 28|. One can argue 
that, 

TcM = rc{<pL = 0) + rcfl(0L), (1) 

where Tc{(f)L) and Tc{(f)L = 0) are characteristic timescales for the motion of a CP in a blend 
and in a pure melt (no LPs), respectively, and tcr is the characteristic CR timescale. When 
(pL ~ 0, it follows that tcr ~ 0, and hence tc{4>l) ~ 'Tcr{^)- Similarly, when Zq ^ 1, we 
expect TcR{<t>L) ^ Tc{(t>L = 0), and tc^cPl) ~ TcBi4>L)- The characteristic timescales for CPs 
and LPs may be approximated by rc(0L) ~ Rc{(j)L)/Dc{4>L), andrL(0L) = Rl{(j)L)/ Dl{4>l), 
respectively. Using standard CR arguments, the local hopping time for the Zc "effective 
Rouse" beads (entanglements on the CP primitive path) is set by the LP relaxation time tl, 
and one obtains the familiar tcr ~ '^lZ^- However, this does not complete the description 
of the CR process, because it assumes that the frictional drag per entanglement segment of 
the CP, or the effective Rouse bead, C,Rouse is a constant. As mentioned earlier, for a CP of a 
given length A^, the number of Rouse beads Zc varies with the composition of the CLE. If C, 
is the monomeric friction coefficient, then the total frictional drag of the CP is C,N , which is 
distributed among the Zc effective Rouse beads. Therefore, Cnouse = C^/^Cj and ignoring 
constants, 

tcr^tlZc ^ tlZc [2) 

Thus, equation [1] can be rewritten as. 



RU^l) RU^l = 0) , ^ Rl 



Dc{<Pl) Dc{<l>L = 0)^^^DL{<l>Ly 
Further, Rl{4>l)/ Rc{4>l) is a weak function of both and (pi, and is easily dominated by 



the change in self-diffusivity [20|, |2l||. Neglecting this change in the size, we rearrange eqn. 
|3]as, 

Dci<PL = 0) _ Dci<PL = 0) . . 

where Ci and C2 are constants that account for prefactors, ignored in this minimal model. 



From our simulations (fig. [T]), and from prior primitive path analysis, 2l| all the parame- 



ters in eqn. Hlcan be determined, and the viability of the minimal model can be ascertained. 
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FIG. 2: To validate equation [H the self-diffusivity simulation data (solid symbols) on = 150 
(circles) and A*" = 300 (diamonds) and experimental data (open symbols) on the self-diffusivity of 
CPs in an entangled linear DNA matrix (ref 25|) at different concentrations, and three different 
molecular weights are replotted. The straight line confirms that there is a linear relationship 
between the quantities on the horizontal and vertical axes. 



As shown in figure [2], all the available simulation data, independent of composition and 
molecular weight, collapse on to a linear master curve. 

To test whether experimental data may also abide by eqn. HJ we employed recently 
published self-diffusivity data on solutions of linear and cyclic DNA [25] . In this study, the 
authors used fluorescence microscopy to measure the diffusivities Dij of tracer LPs and CPs, 
in a matrix of either LPs or CPs, where the subscripts i and j represent the topologies of the 
tracer and matrix molecules, respectively. These diffusivities were compiled as a function of 
contour length (5.9-kbp, 11.1-kbp, 25-kbp and 45kbp), and total solution concentration (up 
to 1 mg/ml), and for small values of these two variables the systems were not entangled. 
Using data on Dll (linear tracer in a linear matrix), we demarcated the transition from 



unentangled to entangled dynamics. Thus, from figures 2 and 3 in ref. [25|, we found that 
the 5.9-kbp sample is too short to be entangled, at any concentration. The concentrations 
at which the shift to reptation dynamics is observed for the pure LPs of length 11.1-kbp, 
25-kbp and 45kbp samples, was 0.7 mg/ml, 0.6 mg/ml and 0.4 mg/ml, respectively. For 



the most well-entangled systems, they observed Dec > Die > Dll > Dql, which is in 
accordance with simulation data on = 300. Thus, there are 16 data-points available in 
the entangled regime at different lengths and concentrations for a tracer CP in a LP matrix. 
Under these conditions, Dc{(t)L) = -Dcl, and 1. Similarly, Dc{(pL = 0) = Dec, and 

DM = Dll. 

The total concentration c, and contour length I contribute differently to the overall dy- 



namics 29|. In the present analysis, the relevant relationships are the dependence of the 
average number of entanglements per chain on these two parameters. We recall that at a 
given c, that the average number of entanglements per polymer chain is proportional to its 
contour length, ~ /. Similarly, for a given /, the entanglement density increases with c. 
Since the modulus G ~ c^^'^, and the entanglement length /g ~ cRT/G ~ c~^/^, it follows 



that Zl ~ l/le ~ C^^"^ [30]. Thus, in the present case, Zc = = {1/Iq){c/ CqY/^ , where 
/o=3 kbp, and Co=l mg/ml were chosen, [25| although it should be pointed out that different 



values for the entanglement molecular weight of DNA solutions at cq =1 mg/ml between 



Iq = 1-30 kbp are supported in the literature |3l|, l32|, |33j . If we superpose these data-points 
(fig. [2]), we find that although two of the 16 data-points diverge from the trendline, the rest 
of the data are in very good agreement with the predicted linear dependence. 

SUMMARY 

We performed Monte Carlo simulations to estimate the self-diffusivity of entangled CPs 
and LPs in blends, and constructed a simple constraint release model. Both, the simulation, 
and experimental data on entangled DNA solutions appear to obey the minimal model over 
a wide range compositions, concentrations, and molecular weights. 
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